Disorder Induced Phase Transition for Random 2D Dirac Fermions 
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We present evidence that two dimensional Dirac fermions in the presence of random Abelian gauge 
potential exhibit a phase transition when the disorder strength exceeds a certain critical value. We 
argue that this phase transition has novel properties unique to disordered systems. It resembles in 
many ways the transition from dilute to dense polymers phase in two dimension. In particular, we 
argue that the central charge of the disordered Dirac fermions, being c = before the transition, 
changes to c = —2. We discuss possible implications to quantum Hall transitions, in view of recently 
proposed model for quantum Hall transitions with c = —2. 
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The study of two dimensional Dirac fermions in the 
context of quantum Hall effect was initiated in Ref. |Q . 
These fermions exhibit an integer quantum Hall transi- 
tion as the value of their mass is tuned through zero. Un- 
fortunately, for pure fermions this transition lies in a uni- 
versality class different from the one observed in quantum 
Hall effect experiments, with the exponents of the Ising 
model. In order to model the more realistic situation, 
various types of disorder should be turned on. Usually 
three types of disorder are studied, random mass, random 
scalar potential and random gauge potential. Random 
mass turns out to be irrelevant in the renormalization 
group sense and the critical properties of the transition 
do not change as we turn it on. Random gauge poten- 
tial results in the critical line, with the critical properties 
of the transition continuously dependent on the disorder 
strength. And the random scalar potential turns out to 
be relevant, flowing away to some unknown critical point. 

It was further argued in Ref. |l| that in order to simu- 
late the realistic quantum Hall transition, all three types 
of randomness should be present. That causes the system 
to flow to a strongly coupled critical point with unknown 
properties. This critical point should be in the same uni- 
versality class as the generic quantum Hall transitions, 
usually described by the Pruisken sigma model (see 
Q for its supersymmetric version). There is very little 
we can say about the properties of that critical point, 
except we might expect the effective field theory which 
describes it to have central charge c = 0. 

In view of the difficulty of approaching the generic situ- 
ation, a number of papers concentrated on the properties 
of critical line generated by the presence of random gauge 
potential jf-||. 

It appears that the correlation functions of such disor- 
dered fermions are very easy to calculate. Indeed, there 
are several methods available, all of them result in the 
correlation functions continuously dependent on the dis- 
order strength 

However, in this letter we will show that previous treat- 
ments were somewhat careless. In fact, the correlation 
functions of the Dirac fermions in the presence of ran- 
dom gauge potential depend continuously on the disorder 
strength up to a certain critical value. If we continue to 



increase the strength of the disorder beyond that value, 
the critical properties of the correlation functions stop 
being dependent on the disorder strength. The system 
as a whole undergoes a phase transition. The density of 
states at zero energy is zero below the transition and is a 
constant above it. And we argue that the effective field 
theory describing that critical line changes considerably. 
Its central charge was c = below the transition, and 
it changes to c = — 2 above the transition. We draw a 
direct parallel with the theory of polymers (self avoiding 
random walks) in two dimensions. They are known to un- 
dergo a phase transition called dilute to dense polymer 
transition. The central charge of the effective field theory 
describing polymers changes from c = to c = —2 [Jij. 
Moreover, we will see that the field theories describing 
the dilute and dense polymer phases and those describ- 
ing the two phases of the fermions with random gauge 
potential are closely related. 

In fact, the first evidence of this phase transition 
showed up in Ref. Q . There the multifractal exponents of 
the zero energy wave function of the Dirac fermions with 
random gauge potential has been calculated exactly. It 
was discovered that these exponents change sharply as 
the disorder strength exceeds a certain value. Evidently 
this behavior is just another manifestation of the phase 
transition discussed in this letter. 

We will conclude with discussing the implications of 
this transition for the not yet established theory of the 
integer quantum Hall transitions. It appears that a natu- 
ral way to approach the generic quantum Hall transition 
point is to turn on strong random gauge potential. That 
generates the constant density of states expected at the 
generic quantum Hall transition, and changes the central 
charge of the theory from c = to c = - 2. After that, 
we may want to turn on random mass and random scalar 
potential, following the logic of jij . The theory will then 
flow to a critical point which, contrary to a naive belief, 
will not be a c = theory. This may shed new light 
on the recent proposal that the physics of quantum Hall 
transitions should be captured by a c = —2 theory ||. 
However, the connection to the model proposed in ^ is 
at this stage purely speculative. 

Now we proceed with the derivation of this phase tran- 
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sition. To begin with, let us fix the notations. We study 
the two dimensional Dirac fermions ip, ip in the presence 
of gauge potential with the Hamiltonian 



H 



d 2 x 4> (id,j, - Afj) cr^ip + eipip, 



(1) 



where a% = a Xl a 2 = <? y are Pauli matrices and e is the 
energy. A^ is supposed to be random with the probability 
density 



P cx exp 



1 



dfx A' 



(2) 



where g is the disorder strength. It is convenient at this 
stage to decompose A^ into a physical part <fi and a pure 
gauge x, following §], 



An, = d^x ■ 



(3) 



where is the antisymmetric tensor. 

In terms of the fields \ an d 'P the probability density 
becomes 



P cx exp 



1 



d 2 x (cU>) 2 + (d„xY 



d 2 xd^0df,9 



(4) 



where the anticommuting fields 8 and 8 provide for the 
Jacobian of the variable change from A^ to <f> and \ ( see 
Ref. @). 

The simplest way to compute the correlation functions 
of (|l|) is to use the bosonization technique. These cor- 
relation functions can be obtained with the help of the 
following path integral 



T>(p exp 



2 x ({d^ipf + ie^A^dytpj 



(5) 



It is very easy to compute the partition function Z by 
completing the square, 



exp 



-\\d 2 x {d^f 



(6) 



which is the effective action of the celebrated Schwinger 
model. 

To compute an averaged correlation function of a phys- 
ical operator, say ip{z) = exp(ip(z)) or any other object 
X, we need to average the following path integral 



From this point on, we could proceed in two different 
ways. One of them would be to calculate the correlation 
function of X and then average it over the random A^ 
(that is, 4>) with the help of (ft). The other way would be 
to average the path integral (ra) over the random A^ first 
and then calculate the correlation function. In a quite 
intriguing way, these two methods do not produce the 
same results. 

Let us first proceed in the former way. We choose 
the left moving Dirac fermion tJj(z) = exp (iip(z)) as the 
object X whose correlation functions we would want to 
compute. We then compute the correlation function of 
exp (itp(z)) by shifting the variable of integration ip by 
the amount proportional to i</>. Then we average over 
the random (f> with the help of (^) (averaging over \ an d 
0, 8 is irrelevant as these fields completely drop out of any 
physical correlation function). As a result, we discover 
that the dimension of the field exp (iip(z)), being equal to 
1/2 without the disorder, becomes 1/2 — 5/8 in presence 
of disorder. Thus we reproduced the standard result of 

However we should become suspicious of this result. 
When g > 4, the dimension of the operator exp (i(p(z)) 
becomes zero and then negative. That means this field 
acquires an expectation value (which, in case of the nega- 
tive dimension, scales with the system size towards bigger 
values). Once it acquires an expectation value, we can no 
longer argue that we could perform the shift in the path 
integral to arrive at this result in the first place. The 
path integral is no longer invariant under the shift of ip. 

We would have arrived at the same problem had we 
analyzed the conformal field theory approach of Ref. || . 
There it was shown how to calculate the dimension of 
ip(z) with the help of U(l|l) Kac-Moody algebra hid- 
den in the supersymmetric approach to (El). In the very 
same way, its dimension becomes negative at large disor- 
der and it acquires an expectation value. However, the 
expectation value for that operator would break the very 
Kac-Moody symmetry which the approach of ]^] utilized 
so successfully. 

A way out of this quagmire would be to add the sym- 
metry breaking energy term of (|l|) to the path integral, 
eipip oc e cos (ip) and try to take e to zero. However, 
adding this term to the path integral (||) would make 
it completely intractable. 

Instead, we should try a different approach. Let us first 
average (||) with the help of the probability distribution 
in other words do the 4> integral in 



(X) = 



JVipXexp - J d 2 x ( (d^p) 2 + ie^A^ip 



Z 



(7) 



over the random A^. With the help of (0) it becomes 



[X) 



T>p X exp 



d x d^ip + -d n 



(8) 



T>(\fDp exp 



d 2 x 



\d^) 2 



(9) 



It is not hard to do that if g < 4 by completing the 
square. That gives back the dimension 1/2 — g/8 to the 
operator exp(iip(z)). 

However, if g > 4, the integral over <p in (JTJ) becomes 
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divergent. A different technique is needed to compute it 
in the strong disorder regime. 

As always when dealing with divergent integrals, it is 
a good idea to regularize them. We are going to limit the 
<p integration by a certain cutoff in the functional space. 

In fact, it is easiest to do that in the ordinary inte- 
gral equivalent to @. Consider the following ordinary 
integral 



dxdy 



■ exp 



— oo J — oc 



2 
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(10) 



It is just an ordinary integral, and yet it captures the 
essential properties of (^). Let us try to compute ( x 2 ) 
with the help of this integral. By that we mean inserting 
x 2 inside the integral in dlO) and then doing the integral. 
If we try to do the integral over x first, and then integrate 
over y, we arrive at ( x 2 ) = 1/2 — g/8. However, if we 
reverse the order, the integral over y becomes divergent 
if g > 4, just like in (Q). So we limit the integral over y 
by the cutoff — L and L. It is not hard to estimate this 
integral at large L, 



dy exp 



1 1 



1 



■ exp 



4 9, 



tyx oc 
cos (xL) 



(11) 



We can now use that estimate to compute ( x 2 ) as in 



dx x 2 exp 



ay l 2 



oc exp 



L 2 
9 



(12) 



It is safe to extend the integration over x to from — oo to 
+oo since it converges very fast. Now it is obvious that 
taking the limit L — > oo we obtain ( x 2 ) =0 for g > 4 
with the term exp (— L 2 /<?) suppressing all others. 

Let us pause now and translate this result to the lan- 
guage of (|9|). It implies that doing the integration over <p 
first, with a suitable cutoff, we obtain that the dimension 
of the field exp (itp) is equal to for all g > 4. 

Thus we arrive at the result announced in the begin- 
ning. The dimension of the field tp of (0) is 1/2 — g/8 if 
the disorder strength g is less than 4, and is equal to 
for g > 4. The field tptp oc cos (y), which is the density of 
states of the fermions, acquires an expectation value at 
5>4. 

Figuratively speaking, integrating over disorder <p sup- 
presses the quantum fluctuations of tp when we are in 
the strong disorder phase. That is why exp (iip) becomes 
zero dimensional. We obtained this result using a clever 
'brute force' averaging of the correlation functions of ([!]). 
However, we should expect to see the same phenomenon 
in the supersymmetry or replica technique. 



In principle we can now use the procedure outlined 
above to compute any correlation function of the theory 
in the strong disorder phase. Therefore for all practi- 
cal purposes we now understand the nature of this phase 
transition. Let us repeat that to see it, we had to impose 
cutoff on the integration over disorder in the functional 
space. The nonanalytic behavior, typical of phase transi- 
tions, has been recovered as we took that cutoff to infin- 
ity. Later we will see that in this respect the transition 
we have been studying has a lot in common with a more 
established dilute to dense polymer transition. 

Now we need to understand what this phase transi- 
tion implies for the effective field theory of the disordered 
Dirac fermions. Effective field theory can be derived in 
either replica or supersymmetry approach, and we are 
going to concentrate on the latter. Supersymmetrization 
of (nj) can be achieved by adding Dirac bosons to it, as 
in 



H = J d 2 xtP (idp - Afj) a^tp + b (id^ - A^) cr M c, (13) 

where b and c are commuting spinors. In calling them b 
and c we follow the accepted convention of ||. 

At this stage we can directly average over A M with the 
help of (||) and obtain the effective supersymmetric field 
theory. 

However, to avoid dealing with a complicated inter- 
acting field theory, it is convenient to bosonize the fields 
tp, tp, and 6, c before the averaging. That is achieved 
with the help of two bosonic fields, <p\ and <p 2 , and two 
fermionic fields 8, 8 as in follows 



(14) 



- f d 2 x ((d„ Vl ) 2 +(d l , V 2) 2 +e l _ ll ,A l _ l (id vVl +a v ip2)-a l ,ed ll s) 



We recall that the field ipi bosonizes the fermions tp and 
tp while the fields tp 2 , 8 and 8 are needed to bosonize the 
b, c fields. 

Averaging over A^ with the help of (||) we obtain 
(dropping the 8 and 8 terms for the time being) 



T>tpiV(f2 e" 



- / d 2 x ((d^^ + ^^f-iiidvi+d^) 2 ) , 15 ) 



It is not hard to see that this path integral and the 
one of (^) are the equivalent, as far as the averages in- 
volving the field ipi or the field tp of (^) are concerned. 
Just as in (||), the integral over Lp 2 in (|15|) is divergent 
if g > 4. Cutting this integral off and then taking the 
cutoff to infinity, one gets the same phase transition as 
we discussed before. Effectively, the integration over ip 2 
at strong disorder freezes out the fluctuations of the field 
tp\. We recall that ip 2 was the result of bosonization of 
the auxiliary fields b, c. The only remaining fluctuating 
fields are 8 and 8 of (|l|). 

It now becomes obvious what this phase transition im- 
plies for supersymmetry. Of course, after the transition 
the invariance under the supergroup U(l|l), manifest in 
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(|l3|), is no longer there. In fact, a good way of think- 
ing about this phase transition is to recall that accord- 
ing to ||, the action Jl3|), after averaging over A^, has 
a Kac-Moody U(l|l) symmetry. It is generated by the 
currents J, j, rj, and fj with the energy momentum tensor 
quadratic in these currents, T oc (Jj+r]fj—fjr])+(4:—g)JJ. 
After the phase transition the currents J and j, being 
just linear combinations of the fields dipi and dif2, anni- 
hilate any physical state. The energy momentum tensor 
becomes T = 7777 and the central charge is now c = — 2. 
The fields rj and fj are nothing else, but the derivatives of 
the fields 8 and 8 introduced in (14). Therefore, the sym- 
metry U(l j 1) changes to the symmetry PSU(1|1) of the 
c = — 2 theory. It is tempting to add that if we worked 
with both the advanced and retarded sectors, the U(2|2) 
symmetry would change into PSU(2|2), in the spirit of 
the recent proposal of Q . 

We are going to change the subject slightly and discuss 
the connection of this phase transition with the physics 
of two dimensional polymers. The polymer problem is 
usually defined as self avoiding random walks, with the 
Green's function given by 



(16) 



where the sum goes over all the nonselfintersecting tra- 
jectories connecting the points x and x' , L is their length, 
and \i is a parameter usually called frugacity (for a re- 
view, see Ref. |10|). When /1 is small, the polymer 
Green's function (|16[) is short ranged. As we increase 
/i, the Green's function becomes critical, and the the- 
ory approaches the dilute polymer critical point. If we 
increase \jl even further, the sum in ( [l6| ) becomes diver- 
gent. At this point, we may want to put the polymers 
inside a finite box, not allowing the trajectories in Jl6| ) 
to go outside that box. Then the sum in ( |l6| ) becomes 
convergent in that new sense, and the theory enters the 
new phase, the so-called dense phase. In this phase, the 
polymers fill the entire space available to them, hence the 
name of the phase [^3| . 

There is a deep analogy between the dilute to dense 
polymer phase transition and the disordered fermion 
phase transition we discussed in this letter. Indeed, in 
the path integral representation of ( |l6| ) we sum over all 
the trajectories connecting the points x and x' , penaliz- 
ing those which intersect. The transition from a dilute to 
dense phase will manifest itself in the divergence of this 
path integral, and we have to cutoff the integral in the 
space of paths (restrict the paths to lie in a finite box) to 
see the dense phase, just like we did with the disordered 
fermions in this letter. 

It is quite instructive to recall that the polymer prob- 
lem can also be cast in the form of a motion in the pres- 
ence of disorder. It is known to be equivalent to the 
problem of a quantum-mechanical particle moving in a 
purely imaginary random potential pl| . Studied by the 
replicated path integral, it can then be mapped into the 
n — ► limit of the O(n) model, as was first noticed in 0] 



and later successfully utilized to find the critical proper- 
ties of the dilute polymer phase. 

Moreover, in a direct analogy with disordered Dirac 
fermions, before the phase transition the critical proper- 
ties of the dilute polymers can be captured by a version 
of the U(l|l) Kac-Moody algebra @Q. After the tran- 
sition, however, the dense polymer phase is described by 
the same c = — 2 theory 0] as the one proposed in this 
letter to describe the strong disorder phase of the Dirac 
fermions. 

Returning to our main problem, the integer quantum 
Hall transitions, we think it is clear now what the fur- 
ther course of action should be. One should study Dirac 
fermions in the presence of random gauge potential and 
other types of disorder. The strong random gauge poten- 
tial effectively freezes out two bosonic degrees of freedom, 
leaving us with a c = — 2 theory with a constant density 
of states. Whether that will result in tractable theory 
is an open question, however. It is also not at all obvi- 
ous that this theory will flow towards a PSU(2|2) sigma 
model with a Wess-Zumino term, as was proposed in [0. 
However, it will flow to a theory similar at least in spirit 
to that sigma model. What kind of theory it will flow to 
definitely deserves further study. 
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